Mon. Not. R. Astron. Soc. 000, 0-0 (0000) Printed 1 February 2008 (MN style file vl.4) 



Overconstrained dynamics in galaxy redshift surveys 



Mikel Susperregi 

Fisika Teorikoaren eta Zientzi Historiaren Saila, Zientzi Fakultatea, Euskal Herriko Unibertsitatea, PO Box 644t 4S080 Bilbao, Spain 
Email: wtpsuxxm@ehu.es 



Version of 1 February 2008 



o 
o 

(N 
0^ 



> 

m 
O 



o 



ABSTRACT 

The least-action principle (LAP) method is used on four galaxy redshift surveys to 
measure the density parameter fl^ and the matter and galaxy-galaxy power spec- 
tra. The datasets are PSCz, ORS, Mark III and SFL The LAP method is applied 
on the surveys simultaneously, resulting in an overconstrained dynamical system that 
describes the cosmic overdensities and velocity flows. The system is solved by relaxing 
the constraint that each survey imposes upon the cosmic fields. A least-squares opti- 
mization of the errors that arise in the process yields the cosmic fields and the value of 
f2„i that is the best fit to the ensemble of datasets. The analysis has been carried out 
with a high-resolution Gaussian smoothing of 500 km s~^and over a spherical selected 
volume of radius 9,000 km s~^. We have assigned a weight to each survey, depending 
on their density of sampling, and this parameter determines their relative infiuence in 
limiting the domain of the overall solution. The influence of each survey on the flnal 
value of ilni, the cosmographical features of the cosmic flelds and the power spectra 
largely depends on the distribution function of the errors in the relaxation of the con- 
straints. We flnd that PSCz and Mark III arc closer to the flnal solution than ORS 
and SFL The hkclihood analysis yields = 0.37 ± 0.01 to la level. PSCz and SFI 
are the closest to this value, whereas ORS and Mark III predict a somewhat lower 
flm- The model of bias employed is a scale-dependent one, and we retain up to 42 bias 
coefficients bri in the spherical harmonics formalism. The predicted power spectra are 
estimated in the range of wavenumbers G.02/iMpc~^ ^ k ^ 0.49/iMpc~ , and we 
compare these results with measurements recently reported in the literature. 
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1 INTRODUCTION 

Galaxy redshift surveys are elements of prime importance in 
the determination of cosmological parameters. Looking at 
redshift-space distortions we are able to obtain maximum 
likelihood estimates of f2m (Fisher, Scharf & Lahav 1994; 
Heavens & Taylor 1995; Baker et al. 1998; Hamilton 1998 
and references therein; Dekel 1999a, b; Hamilton, Tegmark 
& Padmanabhan 2000) and the galaxy-galaxy power spec- 
trum (Fisher et al. 1993; Feldman, Kaiser & Peacock 1994; 
Lin et al. 1996; Sutherland et al. 1999; Hamilton & Tegmark 
2000). The dynamics of cosmic velocity flows is also another 
way to measure Sim (Willick et al. 1997; Dekel, Burstein & 
White 1997; da Costa et al. 1998; Kashlinsky 1998; Sigad et 
al. 1998; Susperregi 2001, hereafter SOI) and the power spec- 
trum (Strauss & Willick 1995 and references therein; Kolatt 
& Dekel 1997; Willick et al. 1997; Freudling et al. 1999; Ze- 
havi & Dekel 1999). In the larger picture, galaxy redshift 
surveys contribute modestly, alongside a myriad of other 
datasets (CMB, SNla data, weak lensing surveys, etc) to the 



undertaking of estimating the key cosmological parameters 
(fim, fio, ^B, or JIq, h, Us, nt, etc)(for a review of how 
this inventory of over a dozen parameters is currently tack- 
led, see e.g. Turner 1999; Primack 2000); all the datasets, 
including galaxy surveys, can be ideally combined by ex- 
ploiting cosmic complementarity (Eisenstein, Hu & Tegmark 
1998,1999) leading to demonstrably concordant, or at least 
not inconsistent, predictions. 

A number of techniques have been employed to maxi- 
mize the parameter information extracted from surveys. The 
Fisher information matrix approach is one example of this 
(Tegmark 1997; Tegmark, Taylor & Heavens 1997; Gold- 
berg & Strauss 1998; Taylor & Watts 2000). This approach 
is sound and will lead to very accurate results when applied 
to forthcoming surveys such as SDSS and 2dF; the compu- 
tational effort can be successfully minimised in those cases 
via data-compression techniques, e.g. via the Karnhunen- 
Loeve "signal-to-noise" eigenmodes (Vogeley & Szalay 1996; 
Tegmark et al. 1998; Matsubara, Szalay & Landy 2000), 
quadratic compression into high resolution powers (Tegmark 
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1997; Tegmark, Taylor & Heavens 1998; Tegmark et al. 1998; 
Padmanabhan, Tegmark & Hamilton 2000), etc. There is, on 
the other hand, a case to be made for combining different 
galaxy redshift surveys and extracting information jointly 
from them, as an ensemble, prior to comparing them with 
different kinds of datasets. As a step prior to exploiting cos- 
mic complementarity, it can be argued that redshift surveys 
are best taken into account as an ensemble than individu- 
ally, however larger and better sampled is any given dataset 
in comparison to its predecessors. It is to be expected that a 
comparison of surveys leaves considerable room to manoeu- 
vre to best rid of, or, at least compensate for, the errors 
inherent in each dataset, stochasticity, unaccounted for con- 
taminants, etc. 

The goal of this paper is to set out a procedure to ex- 
tract an accurate estimate of Om from galaxy redshift sur- 
veys, when the datasets are studied as an ensemble. The 
measurement of Om will then be the optimal estimate for 
all the surveys considered. The method will also enable us 
to estimate the matter and galaxy-galaxy power spectra. 
The backbone of the procedure is essentially the least-action 
principle (LAP) reconstruction employed in SOI. That pa- 
per showed that the LAP method is efficient to measure flm 
from a given galgixy redshift survey and that it conveniently 
breaks the degeneracy between Qm and the bias, thus per- 
mitting to measure f2m within fairly arbitrary bias schemes. 
In this paper we consider a bias relationship that is scale- 
dependent in order to investigate the dynamics more realisti- 
cally. Bias is in any case stochastic as well £is scale-dependent 
(Pen 1998; Tegmark & Peebles 1998; Dekel & Lahav 1999; 
Tegmark & Bromley 1999; Taruya 2000), and the bias re- 
latioiisiiip is fully determined by the distribution P(5g|(5m) 
(Dekel & Lahav 1999; Sigad, Branchini & Dekel 2000). The 
stochastic ingredient is largely sample-dependent, and the 
local scatter that is the measure of stochasticity in P{5g\5m) 
is a characteristic of each survey, its shot noise and environ- 
mental effects. The whole point of studying an ensemble of 
surveys is to minimise this effect, and therefore for simplic- 
ity we neglect the existence of stochasticity in the bias and 
we focus on its scale dependence. 

The galaxy redshift surveys that we examine in this pa- 
per are PSCz, ORS, Mark HI and SFI. The LAP method 
is initially applied to each survey individually in the man- 
ner of SOI with a scale-dependent bias model. Hence one 
obtains the cosmic fields and a measurement of fim for each 
survey. Each survey is regarded as a constraint on the cos- 
mic fields, and as they are all indeed different, the system 
that results from considering all constraints simultaneously 
is unavoidably overconstrained. Strictly speaking it is also 
incompatible. The method to find an optimal solution that 
is a meaningful representation of the ensemble entails re- 
laxing all the constraints. The errors that this relaxation 
brings in are dealt with the least-square method, and this 
yields the best fit. The procedure gives us also an estimate 
of Om, which is really a joint measurement for all the sur- 
veys in the ensemble. From the cosmic fields and the bias 
relationship it is then straightforward to compute the power 
spectra, higher-order correlation functions, etc. 

The article is structured as follows. A brief description 
of the LAP method is given in §2 and the algebra is relegated 
to the appendices, though the necessary steps to undertake 
the numerical implementation are set out in §2.1; the char- 



acteristics of the surveys that are relevant to this paper are 
summarized in §3; the procedure for finding LAP solutions 
for the ensemble of datasets is described in §4 for the most 
general case; in §5 we discuss the application of the method 
on the four surveys in question and their results; finally §6 
discusses the main results of this paper, in comparison with 
similar measurements in the literature. 



2 LAP PRELIMINARIES 

The three ingredients of the LAP method are: (a) two 
boundary conditions and (b) the Lagrangian of the self- 
gravitating matter field. The two boundary conditions 
largely determine the cosmological model, as they are the 
endpoints in the evolution of our Universe. Initial homo- 
geneity, as measured by the Sachs- Wolfe constraint on the 
CMB, suggests the first is 



<5(t^0,x) wO, 



(1) 



where 5 is the matter density constrast. The statistics of the 
primordial S is thus unconstrained. The second boundary 
condition is the Universe as is represented in the current 
epoch by a galaxy redshift survey; this is characterised by 
the galaxy number-count n(s) (s denotes redshift coordi- 
nates, (cz,6,ip)). This constraint is given by (SOI) 



Ps{s) 



gals 



V 



l+ffo(x) 



l-Hao(x)" 



(2) 



where ps relates to the galaxy number-count via ps = 
dn(s)/dsdJ7 (where dO denotes a unit of solid angle), V 

is the volume of the survey, x is the radial comoving dis- 
tance, a is the velocity potential (i.e. v = Va), the prime 
denotes the line-of-sight derivative d/dx, g = (5gais is the 
galaxy number density contrast and the subindex denotes 
the present time. We use comoving quantities consistently 
throughout the paper. The bias model that we shall adopt 
is scale-dependent, such that 



go{k)=b(k)S„{k). 



(3) 



The method is used following the same procedure as that ex- 
pounded in SOI, with the addition of a scale-dependent bias. 
The revelant LAP equations are summarized in Appendix 
A for reference. 



2.1 Numerical resolution 

We apply the algorithm as follows. 

• The dataset "V is transformed into a z-space field Ps(s), 
by computing the discrete derivatives of the galaxy number- 
count n(s), binned to the smallest redshift spacing Az, and 
then the result is smoothed choosing an appropriate smooth- 
ing length Vs ■ We shall only implement Gaussian smoothing, 
W{k) = exp(-fcV?/2). 

• The goal is to find a best fit for the modes 5y{t), ay{t) 
(where y = rim in the spherical harmonics and Bessel func- 
tions decomposition, following the formalism of Appendix 
A). A consistent method to achieve this is to get progres- 
sively closer to the real solution by trying Ansatze of increas- 
ing non-linearity in successive iterations. The starting point 
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is the linear solution, which is approximated by inverting 
the relation 5s oc — V^ao, where 



5s = (Smax/47rA'"gals)ps — 1- 

For the linear 5 field, we estimate 
5o cx (5s(x + xaj)). 



(4) 



(5) 



In the surveys where the radial velocities are part of the 
observational input, in our case Mark III and SFI, the Vr 
data are used to construct the first Ansatz for Qq, which is 
then used to compute 

• The linear fields i?n.o,Q y.o a re the first Ansatz that we 
input in the LAP system ( |A1E| ). We transform these first 



into the coefficients 5y"\ ctj/"', using (AfS), (Af3). These are 



the coefficients of the polynomials of order that approx- 
imate the fields. This initial Ansatz is equivalent to linear 
theory and it permits us to solve the homogeneous system. 
It requires an initial guess of the value of flm , as the matrix 
coefficients depend on it. It is worthwhile to experiment 
with a ample range of values such as O.I ^ flm ^ 1.0, to ex- 
amine the impact of these variations in the convergence of 
the solutions. The solution obtained from the homogeneous 
system is least-square fitted to (jf|) , (^) . This in turn requires 
an bias model, and from (^) using Parseval's theorem (i.e. 
(/•(?) = (/) • (<;) , where (/) denotes the Fourier transform of 
/), we have that (?(x) = b{x)5{x). In the spherical harmonics 
and Bessel function decomposition this means 



9{t, 



y 



hrl5y{t)jl(krX) Ylm{9, if) ■ 



(6) 



This assumes that the effect of bias reflects on a superposi- 
tion of radial modes that are rotationally invariant. In this 
paper we have fitted bri for r < 5 and / < 6. Our starting 
point Ansatz in (^) is hri — 1.0. 

• Th e resulting Sy,ay are then used to construct the RHS 
of ( A19 ), the inhomogeneous system. These quadratic terms 
represent the coupling among the normal modes of the per- 
turbations, and their amplitude is a measure of the degree 
of non-linearity present. 



• Successive iterations thus solve (Af9) and we com 
bine this operation with the least-square fitting the solu 
tion to the constraints ([l|),(^). The procedure eventuall; 
yields the correct Sy, ay if the errors in satisfying 
are monotonously decreasing after each iteration. We use 
variations around the solution to construct the Fisher infor- 
mation matrix 

92£(5,a) 



89, 89, 



= - dp 



3^ 



lnps{5,a), 



(7) 



where = ^^m and 9i = bri for i > 0, conveniently ordering 
the bri in a sequence. The likelihood function L, is defined 
as L = exp(— £). (Q) enables us to compute the relative 
likelihood of estimates of the parameters, particularly Qui- 
The working assumption is that 5s in is a random field, 
resulting from the galaxy number-count, which is modelled 
as a 3D Poisson distribution. This enables us to compute £. 
The strength of the LAP method is that it does not require 
a prior, such as the CDM power spectrum, to estimate the 
likelihood of the parameters. 



3 DATASETS 
3.1 PSC^ 

The sample contains 14,819 galaxies, within a spherical re- 
gion of radius Xmax <^ 45,000 kms~^(its median is at ~ 
8, 500 km s^^and the galaxy number-count tails off sparsely 
up to ^ 30,000 kms~^), extracted form the IRAS Point 
Source Catalogue (Saunders et al. 2000), following the se- 
lection procedure given by conditions log(/60//25) > 0.5 
and log(/I00//60) < 0.75 (where fa denotes the measured 
flux at a microns). The former discards all stars whilst the 
latter excludes galactic cirruses located away from the Zone 
of Avoidance (ZoA). The result retains most galaxies plus 
about 300 mixed contaminants, which are excluded on a one- 
to-one basis. The source catalogue contains 17,060 sources, 
of which 1,593 are discarded as they are either within the 
Milky Way or repeated entries; a further 648 are rejected 
as very faint galaxies and unidentifiable sources. The re- 
sulting 6,500 are a combination of IRAS 1.2 Jy and QDOT 
galaxies plus 3,000 additional sources. The sky coverage is 
~ 84% (the ZoA spans a non-axially symmetric region of 
^ 5°). We adopt the unmasked survey, therefore max- 
imising the sky coverage and a selected spherical subvolume 
of radius a;max ~ 15,000 km s~^. General descriptive infor- 
mation and access to data are available at the PSC z webpage 
http: / /www-astro, physics. ox. ac.uk/~wjs/pscz.htmf 



3.2 ORS 

The sample contains 8, 457 optically selected galaxies, of an 
angular size ^ 1.9 arcmin, located within an approximately 
spherical region ~ 18, 000 km s"^ (Santiago et al. 1995; Hud- 
son et al. 1995; Baker et al. 1998). The dataset is densely 
sampled within a;max ~ 8,000 kms~^, which is the spher- 
ical subvolume that we shall consider, and it is sparse in 
the region x ^ 8,000 kms~^ The sample is a compila- 
tion of optically-selected galaxies extracted from the Upsala 
General Catalogue (UGC; 5 > —2.5°), European Southern 
Observatory (ESO; 5 < —17.5°) catalogue and Extension to 
the Southern Observatory Catalogue (ESGC; covering the 
remaining region just south of the Celestial Equator. The 
sky coverage of these samples is respectively 34%, 19% and 
9%. ORS shares a number of galaxies in common with PSCz 
though it contains a larger fraction of spirals. The predicted 
velocity fields in ORS and PSCz are also in good agreement 
(Baker et al. 1998)). The sky coverage is 62%, i.e. signifi- 
cantly more limited than in PSCz (and even IRAS 1.2 Jy), 
due to the extinction in the ZoA, and it spans over \b\ <^ 20°. 
In the sample utilized ORS is filled with 823 IRAS galaxies 
in order to achieve |fe| ^ 5° coverage. A drawback of ORS 
is its poor uniformity, as extinction affects optically-selected 
galaxies more than it does IRAS galaxies. Further descrip- 
tive information and access to data are available at the URL: 



http: / /www. astro, princeton.edu/~strauss/ors/ 



3.3 Mark III 

The sample contains radial velocities of ^ 3, 400 galaxies 
(Willick et al. 1995,1996,1997), compiled from several sets of 
SO spirals and ellipticals. The sky coverage is the entire sky 
with the exception of the ZoA, at |6| 20° - 30°. The data 
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are located within a spherical region of radius a;max ~ 6, 000 
km s^^that is well sampled, though it is fairly anisotropic, 
as it spans to Xmax ~ 8, 000 km s~^ in some directions and 
is consigned to a;niin ~ 4, 000 km s~^ in others. We consider 
the spherical subvolume of radius Xmax8,000 kms-\ Dis- 
tances are inferred via TuUy-Fisher and _D„ — a distance 
indicators and they entail an error 17 — 21%. The data must 
be carefully prepared to correct for Malmquist biases (fol- 
lowing the procedure set out in Sigad et al. (1998) and used 
by Zaroubi, Hoffman & Dekel (1999); also in a improved 
version by Dekel et al. (1999)). This results in the correc- 
tion of the distances of 1,200 objects. Further descriptive 
information and access to data are available at the URL: 



http://redshift.stanford.edu/Marklll/ 



Table 1. Sampling factor in the surveys 





3? max 


-^gals 


V 


PSCz 


15.0 X 10^ 


14.8 X 10^ 


0.11 


ORS 


8.0 X 10^ 


8.3 X 10^ 


0.41 


Mark III 


6.0 X 10^ 


3.4 X 10^ 


0.39 


SFI 


7.0 X 10=* 


1.3 X 10=* 


0.09 



aimax is the radius of the selected spherical volume in each survey 
in units of km , A^gais is the number of galaxies and -q is defined 



by Q. 



The solution that satisfies (^) and the LAP equations ( ^19) 
is the best fit to the surveys that one can find. 



3.4 SFI 



4.1 Domain of application 



The sample contains radial velocities of ~ 1,300 Sbc-Sc 
galaxies, with inclination ^ 45° north of (5 < —45° and 
galactic latitude |6| ^ 10° (da Costa et al. 1996,1998; Haynes 
et al. 1999a, b). The data are distributed within an approx- 
imately spherical region of radius a;max ~ 7,000 kms~^. 
Individual distances are computed in a similar manner as 
in Mark III, and the distance errors are 15 — 20%. A com- 
parison with Mark III shows a good agreement (da Costa et 
al. 1996; ditto for the POTENT reconstruction analysis of 
Dekel et al. 1999) and it is also in agreement with the pre- 
dicted velocity field of IRAS 1.2 Jy (da Costa et al. 1998). 



4 METHOD 

The LAP method applied to one survey is completely deter- 
mined by the following two constraints: 



Xo{5) = 0, 



,(5,a) =0, 



(8) 
(9) 



where (|8|) is the condition of homogeneity (j^ , and (^j con- 
tains the data given by the survey (b). Any additional 
dataset adds one constraint of the type (9[). Let us suppose 
we have such datasets. Obviously in practice all samples 
contain systematic and random errors and not all A''-f 1 con- 
traints are satisfied exactly by the same solution; therefore, 
the system Xi = with i — 0,1, . . . , N need not have a con- 
sistent solution 5,a. To find a solution, the overconstrained 
system must be relaxed so that the solution satisfies the con- 
straints (^ within the least margin of error, ~ ^i, (* > 0), 
and we optimize this by computing the stationary quantity 



N 



= 0, 



(10) 



where 5 denotes a variation, not a density contrast, and Wi is 
the weight given to each survey in the sum. From the central 
limit theorem, we assume that ti are normally distributed 
about zero, i.e. (ti) — and B(ef) ~ Ci. We will there- 
fore adopt the maximum likelihood values Wi = cr~^, where 
a} is the variance of the distribution P[xi]- The likelihood 
function is then 



L = (27r) 



-N/2 



IV 

in- 



exp 



N r, 



(11) 



The first consideration is to determine the radius of the 
spherical volume to which the LAP solutions will be con- 
signed. Ideally, one would wish to use the largest selected 
subvolume possible, in this case that of PSC2;, ajmax ~ 
15, 000 km s~^. This would nonetheless give too little weight 
to the other datasets, as the constraints Xi ~ given by the 
smaller surveys would have little impact in the outcome of 
S,a (the fractions of the selected volumes of ORS, Mark III 
and SFI, with respect to PSCz are, respectively, 0.15, 0.06 
and 0.10). On the other hand, x^iax should not be too small 
either, as a large part of data from the larger samples may 
thus be left unused. An criterion to determine x]^^ is the 
following. We assume that the weight each dataset ought to 
have in determining the volume of the LAP region is given 
by the density of its sampling, given by the ratio 



(12) 



This measure of sampling is given in Table 1 for the surveys 
of our study. It shows that ORS and Mark III are the most 
densely sampled datasets and SFI is the sparsest. PSCz is 
also on the whole fairly sparsely sampled. The volume of the 
LAP region is hence the weighed sum 



(13) 



Using the values given in Table 1, for the case at hand this 
results in 8, 860 km s~^ . In fact this seems to be a good figure 
also for PSCz as the median redshift of that catalogue is at 



8,500kms-\ 
km s~^. 



We shall adopt the round figure 



9,000 



4.2 Solutions 

We proceed as follows. First, we apply the LAP method on 
each sample separately on their own domain of application, 
following the numerical procedure described in §2.1. Hence, 
for each survey we obtain Si ,ai within a spherical volume of 
radius aj^ax (using the values given in Table 1). This also 
yields an estimate of flm and an evaluation of the bias co- 
efficients bri- The results of the measurement of flm in each 
survey are given in the next section. 

Let us consider these solutions, Si,ai, to construct an 
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Table 2. Likelihood of flm 






o 


^PSCz 


"ORS 


r\/J ^ >^1..TTT 

-* iViarkiii 


Pgpi 


0.20 


0.00 


0.00 


0.00 


0.00 


0.22 


0.00 


0.00 


0.00 


0.00 


0.24 


0.00 


0.02 


0.00 


0.00 


0.26 


0.00 


0.31 


0.02 


0.00 


0.28 


0.05 


0.62 


U.oy 


u.uu 


0.30 


0.42 


0.89 


0.70 


0.08 


0.32 


0.73 


1.00 


0.91 


0.38 


0.34 


0.92 


0.93 


1.00 


0.71 


0.36 


1.00 


0.74 


0.93 


0.94 


0.38 


0.91 


0.33 


0.72 


1.00 


0.40 


0.70 


0.06 


0.40 


0.91 



Table 3. Bias parameters b^i 



The value of the likelihood distribution is normalized to unity at 
its maximum in each survey. We have binned the range — 
0.20 — 0.40 in AHm = 0.02 intervals and computed the solutions 
and the likelihood at each step. 

initial Ansatz of the overall solution that satisfies (|l^) . This 
Ansatz is given by 



N 

E 

j=i 



(14) 



(15) 



Thus we adopt the sampling factor as the best weight to 
average the fields. A few trial and error experiments show 
this is a good starting point. In order to consign this Ansatz 
to the domain (p^), if x]^^ < x'^^ then one discards from 
the survey i all the data that lie within the shell ar^ax < 
X < a;maxi ^Lud if x]^^ > sj^ax then ones makes the fields 
identically equal to zero within the range xl^^^ < x < x]^^, 
effectively "padding" that shell to match the LAP volume. 
In the example of this paper, PSCz data are discarded in the 
range 9,000 km s"^< x <15,000 km s"\ and ORS, Mark III 
and SFI are padded with zeros within the shell xl^^x < x < 
9,000 kms"\ 

The Ansatz ([l^),([l5|) is our starting point to compute 
the LAP solution over the domain (jl^) . The numerical prob- 
lem is then to solve ( |A19| ) , introducing variations in the fields 
5, a and in the free parameters fim and bri to eventually pin 
down the stationary solution in (p^). The numerical resolu- 
tion is carried out following the procedure set out in §2.1, 
where all references to the constraints (jl|), 1^) are substi- 
tuted by (0). 



5 RESULTS 

We apply the method described in §4 to the four datasets in 
question, PSCz, ORS, Mark III and SFI. For each survey we 
construct the galaxy number-count density field Ps(s) using 
a high-resolution Gaussian smoothing length of 500 km s~^. 
First of all, the LAP equations are solved for each survey fol- 
lowing §2.1. This results in an estimate of the cosmic fields S, 
a. We carry out a likelihood analysis to compute the value 
of in each case, and the same analysis also yields an 
estimate of the parameters bri ■ In order to do so, the under- 
lying assumption is that the errors in the measurements of 





r = 


r = 1 


r = 2 


r = 3 






I = 


0.98 


1.01 


1.05 


1.09 


1.12 


1.13 


I = 1 


1.02 


1.08 


1.11 


1.12 


1.15 


1.18 


I = 2 


1.09 


1.11 


1.14 


1.17 


1.20 


1.23 


I = 3 


1.16 


1.18 


1.20 


1.23 


1.26 


1.30 


I = 4 


1.22 


1.26 


1.29 


1.31 


1.33 


1.36 


I = 5 


1.30 


1.32 


1.35 


1.39 


1.41 


1.44 


I = 6 


1.38 


1.40 


1.42 


1.45 


1.46 


1.52 



Estimated values of b^i in the range 0<r<5, 0<i<6. The la 
errors in these measurements are Ki 0.01 (in the worst case, that 
of 643, it is 1.2 X 10~^, so this is the upper bound in the error 
bars of these estimates). 



the fields and the parameters, A^^"', Aal,"'' , AS7m and Abri, 
follow a multivariate Gaussian distribution. However, we do 
not assume that the fields themselves are Gaussian (not at 
the present time, though we find that they are Gaussian at 
the initial time, to great accuracy, and therefore it is safe to 
assume Gaussian initial conditions in the likelihood analysis 
though our formalism does not require us to do so). The 
results of P(r2ni) for each survey are given in Table 2. The 
maximum likelihood estimates for Qm are fim^'"'^ = 0.36, 
n°^^ = 0.32, = 0.34, n!f' = 0.38. The la errors 



„(") 



in these measurements are of the order ~ JF." 



respectively. 



-1/2 



which re- 



5.1 Measurement of flm 

The results obtained above are used to construct an initial 
Ansatz (0), (|l5| ) over the domain (|l^), that we use to solve 
( |A19| ) subject to (|l^. We follow the method set out in §4 
and get the LAP solution that is the best fit to the four 
surveys. The likelihood analysis of the solution yields an es- 
timate f2ni = 0.37 ± 0.01 within the la level. This result 
is quite surprising in view of Table 2, due to the fact that 
the two surveys with the largest rj, ORS and Mark III, seem 
to pin down Qui at a lower value, in the range 0.32 — 0.34, 
and yet PSC2 and SFI predict a Qm that is closer to the 
final value. In this case, the optimization of the morpho- 
logical features of the overdensity and velocity fields in the 
LAP solution with those of the four surveys is such that it 
favours a larger Indeed the surveys PSCz and Mark III 
are closer to the final solution than ORS and SFI. The lat- 
ter two yield the highest xf in (0)i though SFI is in better 
agreement with the LAP solution than ORS and therefore 
XsFi is smaller. PSCz and Mark III have the smaller vari- 
ance of the four, and therefore this causes that they have 
a greater weight in determining the cosmography and the 
cosmological parameters in the solution. 

The likelihood analysis for the bias parameters yields 
the estimates shown in Table 3. We have computed 42 co- 
efficients, truncating the indices at r = 5 and 1 — 6. The 
higher / enables us to probe into the nonlinear scales, though 
a stringent limit is set by the smoothing scale of 500 km s^^. 
The high / coefficients are greater than those in the linear 
range (low I), showing greater bias at smaller scales, a trend 
that is also present in all four surveys when the LAP method 
is applied to each one individually. 
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5.2 Matter power spectrum 

Hereafter we use S^^^ and a^^^ evaluated at the maximum 
likelihood values of the parameters, i.e. f2m — 0.37 and brt 
given in Table 3. We compute the Fourier modes 5(k) of 
^LAP^ The matter power spectrum P{k) is then given di- 
rectly by 

(5(k)r(k')) = (2^)'P(fc)5(k~k'), (16) 

assuming that the modes (5(k) are isotropic and homoge- 
neous, and 5 on the RHS is a Dirac delta function. Strictly 
speaking the fields ought be smooth, though we convolute 
the solutions with W{k), in order to discard residual noise 
generated in the derivation of the LAP solution. The results 
for P{k) are shown in Table 4. The galaxy-galaxy power 
spectrum is given by 

P,,{k) = b{kfP{k), (17) 

and the corresponding results are shown in Table 5. We 
have consigned the analysis to the range of wavenumbers 
0.02 /iMpc"^ ^ k ^ 0.49feMpc"^ that we have binned in 
24 equidistant segments in the log scale. The non-linear part 
of the power spectrum is located at fc ^ 0.3 /i Mpc~^, and it 
is in this range that we expect the LAP method to yield sig- 
nificant information. The resulting velocity dispersion is es- 
timated cr„ fa 286±32 km s^^and af = 0.64±0.03. A COBE- 
normalized flat ACDM model, untilted and with Q^n ~ 0.37 
and A = 0.63, would fit the measured Pgg for F = 0.21 ±0.03 
within the la error margin. The errors in Tables 4 and 5 are 
la. In the case of Table 5, they are larger as they result from 
the product of the errors in P{k) and b{k). In addition to 
the la errors that are the result of the likelihood analysis, 
another source of error is due to the limited number of 5y 
modes, which is propagated in the coordinate transforma- 
tions 5y 5k and bri b{k). These amount to ^ 10% of 
the total error, and have been accounted for. 

Other power spectra of interest are (see e.g. Pen 1998; 
Tegmark & Peebles 1998; Dekel & Lahav 1999) the galaxy 
velocity power 

Pg4k) = r{k)b{k)fP{k) (18) 
and the velocity- velocity power 

P„„(fc) = fP{k) (19) 

where \r{k)\ > 1 is a galaxy-velocity correlation coefficient, 
/ ~ is the linear growth rate. Typically r ~ 1 (Hamil- 
ton, Tegmark & Padmanabhan 2000), though it may depart 
from unity in the nonlinear regime. Using the measurement 
of ilm given in §5.1 and Table 4, it is straightforward to 
compute (^). The resulting errors result from the product 
of AP(fc) and AQ„, = 0.01. 



6 DISCUSSION 

The measurement of fim that we have obtained by apply- 
ing the LAP method simultaneously on PSC2:, ORS, Mark 
HI and SFI is = 0.37 ± 0.01. This is consistent with 
the estimates that we have obtained from each survey sep- 
arately (see Table 2), though ORS and Mark HI predict a 
somewhat lower value, Qm ~ 0.32 — 0.34. The main factors 
in the determination of f2m in our analysis are chiefly two: 



(1) the variance of the distribution P[xi] defined in §4, (2) 
the errors x?, that reflect the resemblance of cosmographi- 
cal features in the survey and the LAP soltuion. The ideal 
combination is a small variance and a minimal value of x? , 
indicating a close resemblance between the survey and the 
LAP solution. The surveys that satisfy these requirements 
best are PSCz and Mark HI. The variance of ORS and SFI 
is not much greater than that of PSCz and Mark HI, so 
these surveys are consistent with the final LAP solution by 
an error of amplitude no greater than 6%. 

These results indicate a higher value of flm than re- 
ported in SOI for IRAS 1.2 Jy {Qm = 0.30 with a linear 
bias of & = 1.1). In order to compare the results of this 
paper with other measurements reported in the literature, 
most of which are in the form of the redshift distortions 
parameter (3 « O!^^ /b, we need to have some information 
about the optimal value of the linear bias in those measure- 
ments. Hamilton (1998) surveys measurements up until mid- 
1997, most of which faU within the region (3 « 0.45 - 0.75. 
A higher measurement of /? = 0.89 ± 0.12 is reported by 
Sigad et al. (1998) from Mark HI, and recently a much 
lower one oi j3 = 0.41155' j2 t>y Hamilton, Tegmark & Pad- 
manabhan (2000) from PSCz and a moderately low one of 
fH = 0.5 ±0.1 by Nusser et al. (2000) from a comparison 
between ENEAR peculiar velocities and the PSCz gravity 
field. Our results are consistent with /3 measurements at the 
lower part of the range, as these allow a density parameter 
in the region ~ 0.3 — 0.4 for reasonable values of the lin- 
ear bias 6 « 1.0 — 1.3. Higher values of /3 require a much 
larger f2m that is not consistent with our results. Such is the 
case of Zehavi & Dekel (1999), who predict f2m ^ 0.6 for 
both Mark HI and SFI, and a value of f2m ~ 0.4 that would 
be consistent with our results has a low likelihood in their 
analysis. 

In Table 4 and 5 we have given the values of the 
matter and galaxy-galaxy power spectra respectively. The 
values of Pgg are consistent, within the la errors given, 
with the results of Sutherland et al. (1999) and Hamil- 
ton & Tegmark (2000) for PSCz. The range of wavenum- 
bers considered by Hamilton & Tegmark (2000) is consider- 
able, 0.01 /iMpc"^ ^ k ^ 300/iMpc~\ and they carry out 
a careful analysis of errors by computing the prewhitened 
power spectrum. Using similar results, Hamilton, Tegmark 
& Padmanabhan (2000) fit their measured Pgg to a COBE- 
normalized, untilted fiat ACDM, Q.^ = 0.3 and = 0.7 
(Eisenstein & Hu (1998)), within a range of wavenumbers 
0.01 /iMpc"^ ^ k ^ 1.00/iMpc"\ Our results is in agree- 
ment with theirs, within the more limited range of wavenum- 
bers that we have estimated, namely 0.01/iMpc~^ ^ k ^ 
0.49/iMpc"\ and our analysis favours a slightly higher 
value of f2m that, within the errors, is fairly close to ~ 
0.4. Furthermore, our analysis predicts a value of the veloc- 
ity dispersion (t„ — 286 ± 32 kms~^, slightly lower than 
one of the values adopted by Sutherland et al. (1999), 
cr„ = 300 kms~^. The analysis of peculiar velocities in 
Mark HI and SFI, by Zaroubi et al. (1997) and Freudhng 
et al. (1999) respectively (see also Zehavi & Dekel (1999)), 
yields a high-amplitude estimate of the power spectrum, 
P{k = 0.1hMpc-^)nln^ = (4.5 ± 2.0) x 10^ (/i"^ Mpc)^ 
Our estimates are consistent with this (from Table 4, we 
have P{k = 0.1 /iMpc"^) ^ (4.63±0.36) x 10^ (/i"^ Mpc)^), 
pinning down the values to greater accuracy. It is also the 



© 0000 RAS, MNRAS 000, 0-0 



Overconstrained dynamics in galaxy redshift surveys 7 



Table 4. Matter power spectrum 



k 

(/iMpc-i) 


P{k) 
(/i-^Mpc^) 


AP(fc) 


k 

(/iMpc-i) 


P(k) 


AP(fc) 


0.0182 


9532 


2266 


0.1010 


4630 


363 


0.0210 


10043 


3451 


0.1165 


3981 


298 


0.0242 


11283 


4016 


0.1343 


3374 


274 


0.0279 


12499 


3220 


0.1550 


2935 


260 


0.0322 


13274 


3061 


0.1787 


1879 


184 


0.0372 


13883 


2874 


0.2062 


1140 


165 


0.0429 


14392 


2249 


0.2378 


826 


149 


0.0494 


12529 


1892 


0.2743 


576 


121 


0.0570 


10302 


1220 


0.3164 


488 


82 


0.0658 


8821 


921 


0.3649 


321 


71 


0.0759 


8427 


887 


0.4209 


267 


55 


0.0875 


8033 


408 


0.4855 


210 


41 



Estimated values of P{k) in the interval of wavenumbers 0.0182 feMpc"^ < < 
0.4855/1 Mpc~ ^ . The range has been binned to produce 24 measurements within this 
interval. The errors, AP(fe), are to Icr level in the likelihood estimates. 



Table 5. Real-space power spectrum 



k 


(/i-^Mpc^) 




k 


(/i-^Mpc^) 


AP(fc) 


0.0182 


4232 


2389 


0.1010 


5013 


408 


0.0210 


9721 


3672 


0.1165 


4521 


353 


0.0242 


11396 


4229 


0.1343 


3882 


299 


0.0279 


13102 


3309 


0.1550 


3621 


275 


0.0322 


13855 


3088 


0.1787 


2288 


192 


0.0372 


14935 


2946 


0.2062 


1720 


174 


0.0429 


16223 


2321 


0.2378 


1211 


165 


0.0494 


13725 


1922 


0.2743 


1088 


132 


0.0570 


11263 


1321 


0.3164 


975 


89 


0.0658 


8860 


974 


0.3649 


720 


80 


0.0759 


8751 


925 


0.4209 


612 


62 


0.0875 


8231 


530 


0.4855 


488 


53 



Estimated values of Pgg . The errors are, like in Table 4, la. 



case that in our estimates of Pgg{k) the error bars are sig- 
nificantly smaller than in Hamilton & Tegmark (2000) and 
Sutherland et al. (1999). 
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APPENDIX A: LAP EQUATIONS 

The cosmological perturbations are governed by (SOI) 

c = i(l+^)v2+a{^+v■[(l+^)v]}-,/<5-i^7-l|v,^|^(Al) 

where (f> is the gravitational potential created by the mat- 
ter density contrast S. The variational equations, namely 
5S/S5 = and 5S/Sa — {5 denotes a variation), yield: 



V = Va, 
and 

5 + V ■ [(1 + <5)v] =0, 



(A2) 
(A3) 

(A4) 
(A5) 



The radius of the Universe a{t) is scaled out alongside h 
with our use of units of km/s for the distance scale. We seek 
solutions for 5,a by utilizing the trial fields 



5(t,x) = ^5.,(t)j,(fc.x)y,„(e,vp), 

y 

a{t,x) = ay{t)ji(krx) YimjO, ip), 



(A6) 



(A7) 



where y = rim. By substituting ([A6|),([A7|) into ([A2|),(|A3[) 

we get 



v(f,x) = ^ ^^a'y{t)ji(krx) + A Jy(a) 



Yimi0,ip), (A8) 



0(t,x) 



f^m^ kr ^ Sy(t)jl{krX) Ylr„{9, If), 



(A9) 



where the coefficients a'y and 3y[ot) are calculated in Ap- 
pendix B. The time dependence of the coefficients 5y, ay 
is modelled by a polynomial of degree that is a sum of 
Chebyshev polynomials, i.e. 



n = l 
JV 



a,(f) = ^4"'r„(t). 



(AlO) 



(All) 



The properties of r„ are summarised in A pp endix C. We do 
not consider the 0th order term i n ( | A10|) , (All) as it is a 
constant term. The polynomials (AlO), ( |A11| ) that permit 
us to approximate 5y, ay, in such a way that the weighed 
sum of the errors in the time interval considered is least, are 
given by the coefficients 



4"' = {Sy\T,,){T,,\T,,y\ 

4"' = (Qj,|r„)(r„|r,o~\ 



(A12) 
(A13) 



wher e (|) is the discrete version of the weighed product 
(^), i.e. 



(A14) 



(/Ifl) = ^ W{tm) f(tm) g{tr. 



© 0000 RAS, MNRAS 000, 0-0 



Overconstrained dynamics in galaxy redshift surveys 9 



the interval [0, 1]. The boundary conditions (SOI) 



where w{t) is (C14) and M is the number of tm divisions on 

(A15) 







n s-(n) 



Ps(s) 



gals 



V 



l+5o(x) 



1 + ao(x) 



(A16) 



go is the present galaxy number-count density contras t and 
V is the volume of the survey. The coefficients 5y enter ( A16) 
through the bias relationship (^), that also introduces the 
coefficie nts bj -i. The evaluation of the modes at the present 



time in ( A16 ) is simplified by virtue of T„(to) — 1. Substi- 
tuting the trial fields into equations ([A4|),([A^), we get 

'^{5y-k'^ay)ji{krx)Yim{^) = -y^< a'y5'y,ji{krx)ji,{kr>x) 

,1 yy' 

X A 3y{S)'\ ■ [x A J^-(a)] I YUn) Yi,^, (fi), (A17) 



+ — 



and 



x)ji'{k,.'x) 



(A18) 



X A Jy{a) 



X A Jyi (a) 



Yi^{n)Yi,^,{n). 



The coefficients Jy{S) are given in Appendix B, as in Jy{a), 
with the substitution a —> 5. By integrating out all coordi- 
nates, it is possible, though arduous, to finally arrange all 
the coefficients Sy"\ in the following inhomogeneous 

matrix system for each y 









r (") 1 
"a 




CO: 


s 

Sy 




V 





y y y 



y" 

(m) ^(m ) 
y" 



(A19) 



where Cy, Cy, Sy, Sy, T^yiyii and fy/^// are computed us- 
ing the o rthogonality relation of the Chebyshev polynomi- 
als ( |Clj), an d the standard orthogonality relations for Yim 
and ji,(p6|), (Df2). These terms are estimated numerically. 



The c oefficients Oy, 5y in the column matrix on the LHS of 
( Af9 ) are arranged vertically in order of increasing n, from 
top to bottom, making a total of 2N entries. The subm atri- 
ces C result from the integration of the LIfS of (|A17|) , and 
the subr natri ces S from integrating the LIfS of ( [Af 8| ). The 
RHS of ( |Af E| ) contains cross coefficients coupling different y 
and n. In conclusion, we have a set of constraints for 5y"\ 



(Af5), (Af6) and (Af9)(as we have int egrat ed out all 



coordinates from the dynamical equations, (A19) is really 
an inhomogeneous system of constraints). 



APPENDIX B: RADIAL DERIVATIVES 

The radial derivative of the velocity potential can be written 
as 



-^a(t,x) = a'y{t) jijkrx) Yi„^(6', ip). 



(Bl) 



where, using the equality 
we have 

' + 1) 



Iji-iiu-) - (' + l)ji+i(w) 



a'y = k, 



I 



(21 + 3) 



{21 - 1) 



(B2) 



(B3) 



Similarly 



„ _ ,2/ (/ + 1) (^ + 2) 

{l + lf f 



+ 



{2l + i){2l + l) (2l-l){2l + l) 
I [l-l) 



(B4) 



(21 - f ) {21 - 3) 



The coefficients 3y{a) in (A8) are therefore given by 



a{l,m + 1) . 

■*y\<^) = 2 "ri(m+l) («xi - X2) 

B{1, m — 1) ... - N . 

H 5 ari(m-i) (2x1 +X2) +imay X3, 



where 
a{l, m) — 

P{l,m) = 



/(/ -I- f) -m(m - f) 



l{l -f- 1) - m{m 4- f) 



1/2 



1/2 



(B5) 

(B6) 
(B7) 



APPENDIX C: CHEBYSHEV POLYNOMIALS 

The Chebyshev polynomials of the first kind r„ satisfy the 
linear homogeneous second-order differential equation 



{l-t )y ~ty + n y = 0, 



(CI) 



where the dot denotes d/dt. Following the convention of 
Courant & Hilbert (f989), they are explicitly given by 



T„{t) 



[n/2] 



2 ^ 

m=0 



(C2) 



m!(n — 2m)! 
that can be rewritten as 

T„(cos6') = cos(n6'). (C3) 
The polynomials satisfy the orthogonality relation 



(T„|T„)=y dtm(t)T„(t)r„(t) =<5,„„-(f + 5„o), (C4) 
where 5„m is a Kronecker delta function and 



w{t) 



2\-l/2 



They have the following property 
■f 



Tn(t)=r:(: 



and their recurrence relation is given by 
2r„(t)T„(f) {t)+T„.^{t), 



(C5) 
(C6) 
(C7) 



for 71 > m. In the asymptotic regime t ^ 00, the zeros of r„ 
are 
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An) 



/2m - 1 

cos IT 

\ 2n 



(C8) 



where t^"^ < t^^^ < . . .. (C4)-(C5) are requ ired in the eval- 
uation of the numerical coefficients in (Alt). 



The interval of orthogonality adopted in (C4) is —1 < 
t < 1. In our application it is more convenient to scale t 
down to the interval [0, 1], and we do so by transforming 
the polynomials via the change of coordinate t — 2t ~ 1. 
Therefore 



r„(t) = r„(2f-i) 



[n/2] 



(C9) 



The coefficients a™ of the new polynomials are found recur- 
sively from the old am through the relations 



m + l ' 

for m — n — l,n — 2, . . 

2 ' 

for m = 0, 1, 2, . . . , n; 



,j; J = 0,1,2, 



for j = 0, 1, 2, . 



(CIO) 
(CU) 
(C12) 



, n; and a 



(m) 



am for m = 0, 1, 2, 



Thus, the orthogonality relation 
becomes 

r-l 



t) in the new domain 



(r„ir„) 



where 



ri;(t) 



dtwit)T4t)Tm{t) = S^mjil + 5„o), (C13) 



-1/2 



-1/2 



(C14) 



Throughout our cal culat ions we use T„ an d the cross- 
product defined by ( C13 ), as opposed to (C4), though we 
have omitted the tildes for simplicity. 



APPENDIX D: ORTHOGONALITY 
RELATIONS FOR Ylm AND Jl 

Dl Spherical harmonics 

A spherical harmonic is defined as (following the Condon- 
Shortley convention) 



Yim{e,v^) 



^ — —Pi (cosc^je ( — 1) 



(Dl) 



47r (l^\m\)\ 

where r = m for m > and r = Q otherwise, and the Pl^(x) 
are the associated Legendre functions, which are given by 

d'+" 



j'r(x) = ^(i-x 



2Nm/2 



/ 2 



1 



(D2) 



Our convention follows Courant & Hilbert (1969) and dif- 
fers from Abramowitz and Stegun (1965) by a factor ( — 1)'". 
( |D2| ) can be otherwise rewritten as 

Pr(x) = i(^(^) F(-M + i;i-^;^),(D3) 

where F is the gamma function and F is the hypergeomet- 
ric function. The associated Legendre functions satisfy the 
orthogonality condition 



pr{x) pr [x)Ax = ^J-^ ii±^ Sir. . 



(D4) 



^ 2l + l(l-m)\ 

The spherical harmonic satisfies the following property 

Y{-^{e,^) = {-irYr{e,H>), (D5) 

and its orthonormal relation is given by 

^27r /"TT 

dif I d{cOSe)Ylm{0,lf)Yi':m'{d,'P) = Sll>Smm'- (D6) 

'0 Jo 
Furthermore, it satisfies the equality 

I 

P(C0S7) = ^^ J2 Ylm{e,>fi)Yim{0',^'), 



(D7) 



where the directions {6,^) and {0',ip') are separated by an 
angle 7, and Pi are the Legendre polynomials, defined as 

p^pr 



D2 Spherical Bessel functions 

The spherical Bessel function of the first kind, j„, is defined 



where J„{z) is the Bessel function of the first kind 
(-1)" 



n[{l + ny. \2 



z] 



(D8) 



(D9) 



that is one of two linearly independent solutions of the 
second-order differential equation 



2 " I ' I /- 2 ,2n „ 

Z W +ZW + (z —I )w = U, 



(DIO) 



where the prime denotes d/dz. One can rewrite (D8) in the 
so-called Gegenbauer's generalization 



Pi{cos9) sinfde. 



(Dll) 



The spherical Bessel function satisfies the following orthog- 
onality relation 

"i 



dxX^jl{krX)jl{ksX) = ^(5rs| J,'_,_l (fcr)^ 



+ 



1 - 



1 + 1/2 



J, 



(D12) 



where Sra is a Kronecker delta function and fcr 7^ 0, and also 
a closure relation given by 



dK K' jl (kx) jl (kx') 



2x- 



■ 5{x — x), 



(D13) 



where 5 is in this case a Dirac delta function. 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science I^'lJrpC style file. 
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